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Abstract 

The properties of the giant-dipole resonance (GDR) are calculated as a func- 
tion of excitation energy, angular momentum, and the compound nucleus 
particle decay width in the nuclei ^^'^Sn and ^''^Pb, and are compared with 
recent experimental data. Differences observed in the behavior of the full- 
width-at-half-maximum of the GDR for ^^'^Sn and ^''^Pb are attributed to the 
fact that shell corrections in ^'^^Pb are stronger than in ^^^Sn, and favor the 
spherical shape at low temperatures. The effects shell corrections have on 
both the free energy and the moments of inertia are discussed in detail. At 
high temperature, the FWHM in i20Sn exhibits effects due to the evaporation 

width of the compound nucleus, while these effects are predicted for ^'^^Pb. 
PACS number: 24.30.Cz 
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I. INTRODUCTION 



The study of the properties of the giant-dipole resonance (GDR) at finite intrinsic exci- 
tation energy has been the objective of many experimental programs during the past decade 
(see the reviews in Ref. These experiments yield important information regarding nu- 
clear motion as a function of temperature. In particular, the role played by quantal and 
thermal fluctuations in the damping of the giant vibrations. In this connection, one can 
individualize the following central issues: (1) the temperature dependence of the intrinsic 
width @J^; (2) the time scale for thermal fluctuations testing the validity of either the adi- 
abatic picture or the occurrence of motional narrowing [^H; (3) the existence of a 

limiting temperature for the observation of collective motion in nuclei |1lO|,P|; and (4) the 
influence of the lifetime of the compound nucleus on the observed width of the GDR |Tl|. Of 
particular importance to address these issues is a systematic and comprehensive comparison 
between experiment and theory over a wide range of temperatures for several nuclei. 

One of the principal experimental techniques for observing the GDR in hot nuclei has 
been compound-nuclear reactions induced in heavy-ion collisions For the most part, the 
wide range of experiments performed so far indicate that the full-width-at-half-maximum 
(FWHM) of the GDR strength function increases as a function of temperature as is pre- 
dicted by theories for the GDR in hot nuclei that account for adiabatic, large-amplitude 
thermal fluctuations of the nuclear shape @J^J^. Many of these experiments, however, in- 
volve slightly different compound systems and are often analyzed using different parameters 
- most notably the level-density parameter. In addition, because of the dynamics of heavy- 
ion collisions, the compound system is generally formed at high angular momentum. Indeed, 
those systems corresponding to the highest excitation energy typically have the largest angu- 
lar momentum content. As such, it is difficult to separate the effects due to large- amplitude 
thermal fluctuations of the shape from those due to angular momentum. 

Recently, two experimental methods for studying the effects of excitation energy and 
angular momentum separately on the GDR have been introduced. In experiments involving 
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compound nuclear reactions, large arrays of gamma detectors have been used in order to 
identify GDR photons associated with a system of definite angular momentum. With this 
experimental setup, the GDR may be studied within an angular momentum window that 
is usually of the order 10-15 units of angular momentum wide, and centered between SO- 
SO h |T2|. An alternative technique is to excite a target nucleus by inelastic scattering with 
light particles ||13|, which, because of the light mass of the projectile, yields an excited system 
with a fairly small angular momentum. By comparing data from these experiments with 
theoretical predictions, it is now possible to analyze the GDR in hot nuclei in terms of the 
effects due to thermal fluctuations and angular momentum separately. 



In an earlier letter |TJ], we presented the results of a systematic study of the FWHM for 
the giant-dipole resonance as a function of temperature, angular momentum, and intrinsic 
width for the nuclei ^^°Sn and ^°'^Pb in comparison with recent experimental data from 



inelastic alpha scattering [|T^]. In this work, in addition to providing the details of how 
this study was carried out, we also expand upon that work by providing a prediction for 
the influence of the evaporation particles on the FWHM in ^°'^Pb at finite temperature. 
Because of the systematic analysis over a range of temperatures and the relatively low 
angular momentum of the emitting system, it is possible to draw conclusions regarding 
the separate roles played by shell corrections, angular momentum, and the lifetime of the 
compound nucleus on the observed width of the GDR. 

This work is organized in the following manner. In section II, the formalism for cal- 
culating the effects of thermal fluctuations of the nuclear shape while projecting angular 
momentum is outlined. A model for the GDR utilizing a quantal, rotating harmonic oscil- 
lator is given in Section III. A description of the shell corrections to the free energy and 
moments of inertia is presented in Section IV, while results and conclusions are given in 
Sections V and VI, respectively. 
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II. THERMAL FLUCTUATIONS 



The description of the GDR in hot nuclei begins by noting that at a finite temperature, 
T, large-amphtude thermal fiuctuations of the nuclear shape play an important role in the 
observation of nuclear properties. Under the assumption that the time scale associated with 
thermal fiuctuations is slow compared to the shift in the dipole frequency caused by the 
fiuctuations (adiabatic motion), the GDR cross section consists of a weighted average over 
all shapes and orientations. Projecting angular momentum, J, the GDR cross section is 
evaluated via |]r^|TB 



where E is the photon energy, = jS'^ d (3 smlSj) d'j sin 6 d6d(j)d'ip is the volume element, 

with a denoting the deformation paramters (3 and 7 and the Euler angles 0, 6, and ip, and 
Zj = /I?[a]/X3/2e-^/^. The factor I{(3,-f,e,ilj) is given by 

J(/5, 7, e, ip) = h cos^ ip sin^ 9 + I2 sin^ ip sin^ 6 + I3 cos^ 9, (2) 

where the Ik represent the deformation- dependent principal moments of inertia. The free 
energy is given by 

F{T, a, J) = F(T, a, u^ot = 0) + ( J + 1/2) V2J(/3, 7, 6, ^), (3) 

where F(T,a, Urot = 0) is the free energy evaluated in the cranking approximation with 
rotational frequency, Uj-ot, equal to zero. 

In many previous works @,PJ^, a different procedure involving a fixed rotational fre- 
quency method for projecting angular momentum has been used. In this formalism, Eq. (|I]) 
would be replaced by 

a{E) = Z-J J V[aHa, E)e-^(^'"'-)/^, (4) 
where = J 'D[a]e~^^'^ and the free energy is given by 



4 



F(r, «, uj) = F{T, a, LUrot = 0) - ^ J(/3, 7, ^, V')^;'. (5) 

In this scheme, the rotational frequency is determined such that the average angular mo- 
mentum of the system is given by P^ , P^ 

( J) = J + 1/2 = T-^ In = Z-'uj I V[a]I{P, 7, 9, V^)e-^(^'"'-)/^. (6) 

The primary disadvantages of the fixed rotational frequency approach are that angular 
momentum is projected only on average and that for finite angular momentum the nuclear 
free energy in Eq. exhibits a saddle point beyond which the system is unstable. This is 
illustrated in Fig. 1, where, in the lower panel, the free energy for ^°^Sn is plotted along the 
oblate noncollective and prolate collective axes at a temperature of 2 MeV and a rotational 
frequency of 1.25 MeV, which corresponds to an average angular momentum of approxi- 
mately 55^. The free energies were computed as described in Section IV, and effectively 
consist of only the liquid-drop component. In the upper panel of Fig. 1, the Boltzman weight 
factor exp[— (F — Ff.q)/T], where F^g is the minimum of the free energy below the saddle 
point, is also plotted. From the figure, it is clear that at high temperature and high angular 
momentum, the presence of the saddle point can be a serious drawback, as it is not possible 
to perform the thermal averaging. In addition, an important shape transition occurring at 
high spin, known as the Jacobi transition, which is characterized by the sudden evolution 
from an oblate noncollective shape to a prolate collective shape with large deformation, is 
absent. The formalism of Eq. (|l]) was introduced in Ref. ||15| to account for these deficiencies 
and to permit a description of the GDR at very high spin. 

In Ref. , an additional method, where only the z- component of the angular momen- 



tum is projected is also presented. In this case, Eq. (|I]) is modified to 

where Zj^ = J P[«]/X^/^e^^/^ and the free energy is given by 

F(T, a, J,) = F(T, a, Urot = 0) + ( J.) V2X(/5, 7, 0, ^). (8) 



The principle feature of this projection method is to give a better overall description than 
Eq. for nonscalar observables such as the angular distribution 02 coefficient, which is 
defined by 

a{E,9) = aiE)[l + a2{E)P2{cose)], (9) 

where 6 is the angle between the observed gamma-ray and the polarized spin direction. In 
heavy-ion fusion experiments, J and lies in a plane perpendicular to the beam direction, 
and 9 is measured relative to the incident beam direction. Then, 02 may be written in terms 
spherical tensor components cr^ of the GDR cross section via 



a2{E) ^ 



ao(^)-^(ai(E) + a_i(E)) 



(10) 



a{E) 
with (T = Y.^iCr^. 

Here, we have performed calculations at low spin using all three methods of angular 
momentum projection, and find that for the FWHM all three methods give the same value to 
within a few hundred keV, with Eqs. (^) and (|I|) giving the largest and smallest, respectively. 
At much higher spins, J ~ 50^, however, Eqs. (|I|) and yield very different results because 
of the presence of the saddle-point barrier in the fixed rotational frequency scheme that does 
not account for the Jacobi transition, and limits the effect of thermal fluctuations. These 



issues are discussed in further detail in Ref. 16 . 



III. MODEL FOR THE GDR 

In principle, the most appropriate description of the GDR strength function in a hot, 
rotating nucleus would be obtained by performing random phase approximation (RPA) 
calculations for each deformation and orientation. Because of the large number of points 
required in performing the thermal averaging of Eqs. (|1|), (^), and (|^), however, this proce- 
dure is computationally impractical. Instead, we make use of the fact that RPA calculations 
indicate that the GDR is a strongly collective excitation that is also rather stable with tem- 



perature . As such, for all practical purposes, the GDR may be modeled by a harmonic 



vibration along the three principal nuclear axes with frequencies inversely proportional to 
the radius of each axis |18|. Variations of this approach (with both quantal and classical 
oscillators), have been used in the past ||^,^,^|^|[, and for completeness, we describe in detail 
the model used in this work in the present section. 

A harmonic oscillator description of the GDR may be derived from a many-body nuclear 
Hamiltonian H with a pure harmonic-oscillator single-particle potential and an isovector 
dipole-dipole interaction as the only two-body term . For the general triaxial nucleus, 

we have 



k=l 



21 



(11) 



i=l \i=l 

where is the third component of the isospin, the oscillator frequencies Uk are inversely 
proportional to the radius along the axis k = 1,2, 3, with 



UJq = [UJ1UJ2UJS 



:o,cooto^Y^^ ^ MeV, 



(12) 



{h = 1), and Kk is the dipole-dipole strength, which empirically is of the order SMul/A. 

In Eq. (pUD, it is possible to introduce a canonical transformation in which H is split into 
two parts. The first describing the intrinsic nuclear degrees of freedom, while the second the 
collective GDR mode, which may be written as 



HD = ki: {pI + Eld 



2' 
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(13) 



where dk is the giant-dipole operator and pk is the conjugate momentum. Using the Hill- 
Wheeler convention the GDR resonance energies along the three intrinsic axes are 

m 



R 

Ek = Eq— = Eoexp 




5 ^ / 27rk' 



where Eq ^ 80^-^/^ MeV is the dipole energy for the spherical shape. 

If the intrinsic nuclear frame is rotating with angular velocity u, Eq. 
modified to include the coriolis and centrifugal forces, becoming 



(14) 



must be 



Hd = \ T.iPl + EM) - u ■ {dx p), (15) 

k 

where u may be taken along the 2;-axis in the external, fixed reference frame, and while pro- 
jecting angular momentum is taken to be the saddle-point value uj = (J+l/2)/X(/5, 7, 9, tfj), 



i.e., the frequency that maximizes the exponential factors in the projection integral |T5|,|16 
In terms of creation and annihilation operators a\ and a^, Eq. (|15D may be written as 



i I E 

Hd = Ekiaiak + afeOfc) + 2 ^ (^ijk^^iJ a]al - a^a^ + ajal - ajUk . (16) 

k ijk V J 

Consolidating the notation, we may write Hd as 

Hd = lYl {^jkalcik + A*f,aka] + Bjka]al + Bji^ajak) (17) 

jk 



with 



Ej — Ek , 
Bjk = ieijkiOi — , (18) 



and 



'^jEjEk 



Ej, a j = k; 



At this point, we note that H^ is only a quadratic function of the coordinates, and, 
therefore, it is possible to introduce a canonical transformation 

01 = Y: {Xla\ - Yl^ak) , (20) 

k 

= E (^r4 - yrak) , (21) 

k 

such that the Hamiltonian may be written as 

Hd = \Y.Eu {Olo, + 0,0l) . (22) 
1/ 

The eigenenergies and transformation coefficients X and Y are found from the 6x6 RPA-like 
eigenvalue problem 

]=eJ ]. (23) 

A* B* J J \0 -1/ {Y" 

8 



Note that the eigenvalues Ei, come in plus-minus pairs, and the three principal modes of the 
GDR correspond to the three positive eigenvalues. 

In order to evaluate the GDR photo-absorption cross section, it is necessary to calculate 
matrix elements of dj. In terms of the creation and annihilation operators Oj, and Oi, we 
have 



and hence the matrix element {v\dj\{]) can be written as 



(24) 



(H^.|0) = ^^(xf + (25) 

In addition, the transition matrix elements must be evaluated in the non-rotating laboratory 
frame. This is accomplished by first transforming the fixed laboratory coordinates to the 
frame rotating about the fixed z-axis with rotational frequency a;, and then into the intrinsic 
frame defined by the Euler angles. To do this, it is necessary to evaluate the matrix elements 
of the spherical tensors d^. Here, we write d^ in terms of its spherical components, that is 

df. = Y.9f.jdj, (26) 

j 

where the matrix g^^j is defined by the well known relations 

f c?3, if /i = 0; 

[T^{di±d2), if/i = ±l. 
The matrix elements in the frame rotating about the 2;-axis become 

'M;io) = E(HrfM'|o)<l(i^) 



where D|j^l(f2) is the rotation function for tensors of rank 1. 



The GDR cross section to be used in Eqs. (|I]) is now readily calculable. From Fermi's 
Golden rule, a{a,uj;E) evaluated in the intrinsic frame for a nucleus with A nucleons, Z 
protons, and neutrons is 
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dir'^p^f) 2ZN 

aUo^, u; E) = — \ {v\d^\^) \l^E [5{E - u)) - 5{E + E,{a, u))] . (29) 



3"^c A 



Noting that 



1 pOO 

6{E -E') = — (30) 
2tx J-oo 



we have 



47rV^2ZA^^ /""^ ,,,,,, m2 ^ 



OO r- 

^i{E-E^{a,Lj))t _ ^i{E+E^{a,oj))t 



An^h 2ZN 



3mc A 



/OO 
dtEe'''' 



(0|4(0)|z/),,^(z/|d^(t)|0),,^ - (0|4(t)|z/),,^(z/|rf^(0)|0),J , (31) 

where in the Heisenberg picture (i^(t) = e~*''^*(i^(0)e*^*. The fact that the experimental 
giant-dipole resonance has an intrinsic width, T^, can be accounted for in Eq. ( pi] ) by 
modifying the exponential by e^*^"'"'^/^)* giving 

a,rA(^,uj;E) = -^^Y.\i^Kml.E 

[BW{E,E,{a,iu),T,)-BW{E,-E,ia,iu),T,)], (32) 

where BW{E, E', T) is a Breit-Wigner function 

BW(S,F.r) = ^^^-^i;^. (33) 

In the limit that u = (i.e. E^, = E^), Eq. ([3^ ) is a sum of three normalized Lorentzians 



each with a centroid at E^ = ^E"^ + TI/A and width Fj,, and satisfies 100% of the classical 
sum rule. For finite u, however, Eq. (|32|) is a sum of three Lorentzians with a normalization 
of the order Ey/Ek, and does not necessarily satisfy 100% of the classical sum rule. 

Lastly, in order to evaluate a{a, uj; E) in the non-rotating laboratory frame, it is necessary 
to evaluate the matrix elements (i^|(i|^''|0) in Eq. (pl|). These matrix elements may be related 
to those in the intrinsic frame via Eq. ( |28D by noting that the transformation from the fixed 
frame to the rotating frame is accomplished by a rotation about the 2;-axis by the angle ut. 
That is, 
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d':\t) = e^^'d'^it) = Y.^^^'d,,{t)Df^,{n). (34) 

From Eq. (pT]) , we see that in addition to mixing the strengths of the various components, 
the GDR energies are themselves shifted by —^uj. Therefore, in the laboratory frame, we 
have 



E[BW{E,E,{a,uj) - fiu,r,) -BW{E,~{E,{a,uj) - fiu),r,)]. (35) 



IV. SHELL CORRECTIONS 

Due to the exponential dependence in Eq. (|l[), the most important ingredient for the 
calculation of the GDR strength function is the nuclear free energy. Here, the free ener- 
gies were computed using the standard Nilsson-Strutinsky procedure extended to finite 
temperature [^, namely 

F = Flo + F^ — Fs = Flo + Fshl, (36) 



where Fld is the liquid-drop free energy evaluated with the parameters of Ref . EM , and F, 



and Fs are the Nilsson and Strutinsky components comprising the shell correction, Fshl, 
to the free energy. In this work, the Nilsson parameters were taken from Ref. ||25|. For 



the most part, the shell corrections for ^^°Sn were found to be quite small (a few hundred 
keV at T ~1.25 MeV), and for all practical purposes can be ignored. This is primarily due 
to the fact that the separate proton and neutron contributions are approximately equal in 
magnitude, but opposite in sign, and, hence, essentially cancel. This is in sharp contrast 
to the strong coherence found in ^''^Pb, where, at low temperatures, strong shell corrections 
(~ —14 MeV at T = MeV) are found that favor the spherical shape. 

We have also investigated the influence of the pairing interaction, and have found that 
effects due to pairing are significant only for temperatures below ~ 0.75 MeV, which is a 
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lower temperature than for which experiments have been performed. In addition, Nilsson- 
Strutinsky calculations that include pairing, indicate that, for the most part, the effects on 
the free energy are negligible. This is because ^°^Pb is a doubly closed-shell nucleus with 
pairing gaps equal to zero for the spherical shape, and in ^^°Sn, as was the case for the 
free energy without pairing disscussed above, the separate proton and neutron contributions 
tend to cancel, leading to a free energy whose deformation dependence is essentially that of 
the liquid drop. 

We note that a numerical determination of the effects of thermal fluctuations in Eq. (|l]) 
in general requires an exploration of the five dimensional space spanned by the deformation 
and orientation degrees of freedom, in which a large number of points are required in order 
to assure sufficient accuracy (especially at finite angular momentum). In this regard, a 
Nilsson-Strutinsky calculation for each point may be too time consuming. Therefore, it is 
useful to parameterize the free energy using functions that mimic the behavior of the Nilsson- 
Strutinsky calculation as closely as possible. It has been pointed out |^ that, being a scalar 
quantity, the free energy must be a function of the rotational invariants of the quadrupole 
deformation, that is 

F(r, 13, 7) = Fo(T) + A{T)[3^ - B{T)[3^ cos(37) + C{T)P^ + ... (37) 

Although this Landau parameterization gives a good overall description of the free energy, in 
particular regarding to shape transitions, it may not be adequate for the evaluation of Eq. (|l]) 
because at somewhat larger deformations Eq. ( pT]) deviates from the Nilsson-Strutinsky 
calculation, often giving a much stiffer free energy. This is principally because Eq. (|37D 
attempts to combine both the liquid-drop free energy and shell corrections, Fshl = F^ — Fs, 
into the same parameterization. An alternative approach is to parameterize instead only 
the shell corrections to the free energy using a function of the rotational invariants. 

Exhibited in Fig. 2 (solid points) are shell corrections to the free energy at Urot = 
as a function of temperature for oblate (7 = vr/3), prolate (7 = 0), and triaxial (7 = 
7r/6) deformations for ^°^Pb. The general overall behavior of Fshl is to decrease with 
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temperature, gradually melting {Fshl ~ MeV) for temperatures of the order T = 2.5 MeV, 
and that they tend to oscillate with deformation, but appearing to be damped at larger jS. 
In this light, it is possible to parameterize Fshl with a series of functions that are in 
fact themselves functions of the rotational invariants /5^cos(37), etc... One possible 
parameterization is 

Fshl{P,1,T) = Adi{Bif3)CiT/smh{CiT) 

1=0 
odd 

+ ^ AaiiBiP) cos(37)aT/sinh(aT), (38) 

1=3 

where the ji are spherical Bessel functions. We note that C/T/sinh(CiT) is the expected 
attenuation behavior as a function temperature when the single-particle Hamiltonian is a 
degenerate harmonic oscillator . 



The parameters Ai, Bi, and Ci can be determined in the following manner. First, carry 
out a Nilsson-Strutinsky calculation for oblate, prolate and triaxial shapes up to /? = 1.0, 
and for temperatures between T = 0.25 and 3.0 Mev. Then fit both parameters Ai and Bi 
to the Nilsson-Strutinsky calculation at T = 0.25 MeV for all three shapes simultaneously. 
Note that at /? = 0, the free energy is completely determined by Aq, and as such is not fit 
upon. In addition, note that the 7 = 7r/6 points are dependent only on the even functions 
in Eq. (^). Typically, the number of terms in Eq. ( ^8]) can be truncated to / ~ 5. With 
the parameters Bi determined at T = 0.25 MeV, these same values are then used to fit 
the Ai values at all other temperatures, giving the sequence {A;(Tj)}, which is then fit to 
the function AiCiT/sinh.{CiT). Shown in Fig. 2 with the solid line are the results of the 
parameterization of the shell corrections to the free energy for ^°^Pb, and the associated 
parameters are listed in Table I. Generally speaking, the parameterization of Eq. ([38| ) gives 
a good overall reproduction of the shell corrections, Fshl, that is rather quick to implement 
with Eq. (|I|). 

We note one feature of the parameterization for ^osp]^ jg ^j^^^^ ^j-^g parameterized shell 
corrections tend to "melt" a little too quickly. For example, for T > 1.5 MeV the parame- 
terized shell corrections underestimate the Nilsson-Strutinsky values by a few hundred keV. 
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It is to be noted, however, that at these temperatures, this amounts to a relatively small 
change in the overall deformation dependence of the total free energy, which, in addition to 
be divided by the temperature in the Boltzman factor, is at basically dominated by 

the liquid-drop component. 

In addition to the free energy, shell structure can also modify the moments of inertia. 
Again, we employ the Nilsson-Strutinsky procedure at finite rotational frequency, and obtain 
shell corrections to the rigid-body moments of inertia, namely 

I = Irigid + In — Is = Irigid + IsHL, (39) 

where here the rigid-body values were evaluated with the radius R = 1.274^/'^. Choosing 
the rotational frequency along the z-axis, the leading behavior as a function of rotational 
frequency for each of the free energy components in Eq. ( P^ ) is given by 

7, T, u) = F{P, 7, T, = 0) - i/3 7, T)u^ (40) 

The moments of inertia I3 can then be obtained by performing a quadratic fit to the free 
energy components. 

In a manner similar to the shell corrections to the free energies, the shell corrections to 
the moments of inertia may also be parameterized by series of Bessel functions, i.e., 

even 

/3^^^(/3,7,T) = Y: Aiji{B^(3)ClT/smh{ClT) 

1=0 
odd 

+ J2AijiiB{p) cos(37)C/T/sinh(C/T) 
1=3 

+ aiJi{Ki(3) cos(7 + 27r/3)r]iT/smh{r]iT) , (41) 
i>i 

where the third term in the sum is included because of rotational invariance arguments for 
the moment of inertia . Once the third component of the moment of inertia is determined 
as a function of T, /3, 7, the remaining two components are obtained by the relations 

/i(T,/3,7)=/3(T,/3,7 + 27r/3), 

/2(T,/3,7)=/3(T,/3,7-27r/3). (42) 
14 



The parameters Aj, Bl , C[ , ai, ki and rji were determined in a similar manner as those for 
the free energy in Eq. (pSD- Again, the shell corrections for ^^°Sn were found to be small and 
negligible. For comparison, both the parameterized and Nilsson-Strutinsky shell corrections 
to the moment of inertia for ^osp]^ ^^^^ shown in Fig. 3 as a function of temperature and for 
the deformations 7 = 7r/3, 7r/6, 0, — tt/S, and — 27r/3. The most important feature is the 
strong shell corrections at /5 = that significantly reduced the moment inertia below the 
rigid-body value. 

Of particular importantance for the moments of inertia is the fact that the spin-orbit 
and terms in the Nilsson Hamiltonian lead to moments of inertia that are approximately 



20-30% larger than the corresponding rigid-body values [0. As such, the shell corrections 
to the moments of inertia used here were reduced by 25%, which corresponds to the average 
difference between the rigid-body and Strutinsky moments of inertia. The corresponding 
parameters Aj, Bl , C/, ai, ki and rji are then listed in Table II for ^'^'^Pb. 

Because of the T~^/^ dependence in the "effective" volume element in Eq. (|I]), it might 
be expected that the strong shell corrections to the moment of inertia in ^°^Pb would sig- 
nificantly affect the GDR strength function, as they appear to give a stronger preference to 
the spherical shape. We find, however, that because of the factor in T>[a], the strong 
shell corrections in X favoring the spherical shape have very little effect on the FWHM of 
the GDR at low spin beyond that produced by the free energy. This is exhibited in Fig. 4, 
where an "effective" weight factor (which for the sake of simplicity ignores the sin 37 factor) 
W = /?^/X^/^e~^/^ is plotted for oblate and prolate shapes at T = 1.0 MeV for various com- 
binations of the shell corrections. In the top panel of the figure, the weight factor is plotted 
including shell corrections to the free energy as well as with and without shell corrections 
to the moments of inertia, whereas the corresponding figure without shell corrections to the 
free energy is shown in the bottom part of the figure. In both cases, it is seen that the overall 
behavior of the weight function is governed by the exponential of the free energy, which is 
plotted in the upper right-hand panel. In addition, the ratio Tlo/Xshl is shown in the 
lower left-hand panel, where it is seen that without the /J*^ factor the spherical shape would 
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have approximately 40% more weight when shell corrections to the moments of inertia are 
included. 



V. RESULTS 

In this section we present the results of a systematic comparison between theoretical 
calculations and recent experimental data |]13| as a function of temperature for both ^^°Sn 



and ^"'^Pb. The thermally averaged GDR cross sections were computed using Eqs.(|T|) and 
(^). In keeping with experimental findings [0], the intrinsic dipole widths, T^, were taken 
to be dependent on the centroid energy via T^, = TQ^Eiy/ EqY , where Eq and Tq are 
the centroid and width for spherical shape and 6 ~ 1.8. The parameters Eq and Tq were 
taken from ground-state data and are Eq = 14.99 MeV and Fq = 5.0 MeV for ^^°Sn and 
Eq = 13.65 MeV and Tq = 4.0 MeV for ^'^^Pb, respectively. Finally, in accordance with the 
considerable theoretical evidence presented in Ref. 0, the intrinsic width Fq is taken to be 
independent of temperature throughout this work. 

Shown in Fig. 5 are the results obtained for the FWHM of the GDR strength function for 
both ^^°Sn and ^'^'^Pb as a function of temperature in comparison with recent experimental 
data. The solid line represents the theoretical values obtained with zero angular momentum. 
The dependence of the FWHM for ^^°Sn and ^osp]^ angular momentum at T = 1.6 MeV 
is illustrated in Fig. 6, where it is seen that for J < 25h the FWHM is essentially unchanged 
from the J = Oh value. Given that the largest average angular momentum in the systems 



studied experimentally is of the order 20h [13|, the effects due to angular momentum on the 



data set of interest are then expected to be negligible. 

As is seen from Fig. 5, theory provides an overall account of the experimental findings. 
The dependence of the FWHM on temperature is quite different between ^^°Sn and ^°^Pb. 
The FWHM in 208p|^ appears to be suppressed at lower temperatures relative to ^^°Sn. 
This is due to the rather strong shell corrections in 208p|^ that favor the spherical shape 
at low temperatures. The affect of such strong shell corrections is to limit the influence 
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of thermal fluctuations at low temperatures, thereby reducing the observed width. This is 
also illustrated in Fig. 5, where the dotted line in the panel for ^osp]^ indicates the FWHM 
assuming no shell corrections. We note that the shell correction effect and the angular 
momentum dependence was also observed for ^^°Ce in Ref. |0. The fact that the adiabatic 
model slightly overestimates the FWHM maybe due to: (1) uncertainties in the extracted 
temperature; (2) the shell corrections being more persistent at higher temperatures than 
predicted by Nilsson-Strutinsky calculations; (3) the fact that the experimental strength 
functions were fit to a single Lorentzian, while theoretically they are obtained from the 
superposition of many Lorentzians; and/or (4) the presence of non- adiabatic effects that 
would lead to a motional narrowing of the FWHM [0. In keeping with point (1) above, one 
can mention that the temperatures inferred from experiment are sensitive to the choice of the 
level-density parameter, and, as a consequence, are uncertain at the level of approximately 
0.2 MeV. 

The FWHM shown in Fig. 5 are essentially consistent with the adiabatic picture for the 
GDR in hot nuclei, and do not present any evidence for the phenomenon known as motional 
narrowing [0,^, which tends to lessen the effects of thermal broadening on the resonance, 
and, hence, reduce the FWHM. As is pointed out in Ref. 0, however, because of a lack of 
reliable theoretical estimates for the time scales associated with thermal fluctuations, the 
FWHM is not sufficient in of itself to exclude motional narrowing. This is particularly true 
when the time scales for (3 and 7 degrees of freedom are much faster than those associated 
with the orientation of the system. In this case, both the response function and the angular 
distribution 02 coefficient are needed in order to make a differentiation between the two 
regimes. 

We note some slight discrepancies between the adiabatic model and experiment for ^^*^Sn. 
To begin with, the FWHM at T = 1.24 MeV is significantly lower than the theoretical 
prediction and is difficult to explain within the framework of the model. This datum seems 
to point to the existence of strong shell corrections that quickly disappear at T = 1.5 MeV, 
which is in disagreement with the expectations of the Nilsson-Strutinsky procedure. At 
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higher temperatures, T ^ 2.8 — 3.1 MeV, the experimental FWHM is somewhat larger than 
the theoretical values, and may indicate a systematic trend to be observed at still higher 
temperatures. Shown in Fig. 7 is the FWHM for ^^°Sn at T = 3.12 MeV as a function of 
the intrinsic width Fq. At this temperature, the experimental FWHM is 11.5 ± 1.0 MeV, 
and we may infer from this datum a value of Fq = 7.7l2!i MeV, as indicated by the solid 
square (11.5 MeV) and open circles (±1 MeV) in Fig. 7. We note, however, that this 
is consistent with the concept that the width observed for the GDR in hot nuclei should 



be increased because of the evaporation of particles from the compound nucleus ||Tl]|. At 
higher excitation energies, the decay rate for particle evaporation increases, and, because 
of the uncertainty principle, the energy of an emitted GDR photon cannot be known with 
a precision better than Ten = T^eJ°^^ + F"^**^^, where Y^'J^-I'"'^^^"-^^'^^^ is the width for particle 
evaporation before and after the emission of the GDR photon. To account for this effect 
in our calculations, we note that the FWHM shown in Fig. 5 are obtained from the full 
response function, which also includes splittings due to the superposition of the various 
intrinsic modes. On the other hand, F^n represents an uncertainty in the GDR photon 
energy due to the lifetime of the initial and final states. As such, Fc„ should be folded into 
the GDR response function by increasing the intrinsic widths via F'^, Fj, + Fc„. In order 
to estimate Fen for ^^'^Sn we refer to Fig. 2 of ref. |]lOl, where Tev is plotted as a function of 
excitation energy for various values of the level-density parameter a (which is conventionally 
defined as a = A/k, MeV~^, and values of Fe„ are plotted for k = 8, 10, and 12). We 
note that at a given excitation energy, T^y exhibits a strong dependence on a. Indeed, at 
= 150 MeV, there is a nearly a factor of three difference between the results for k, = 8 
and 12. This rather strong dependence on a is considerably diminished, however, when 
converting to temperature, defined as Er^ — Egdr = aT"^, as is shown in Table HI where 
Ybefore{after) are given as a function of temperature for k = 10 and 12. Only at 

the highest temperatures (~ 3.5 MeV), where k, is expected to be of the order 12-13, is the 
difference much greater than a few hundred keV. Taking k = 12, we deduce at T ^ 3.1 MeV 
Fen ~ 2.1 MeV, which is in good agreement with the experimental results as is illustrated in 



Fig. 7. To further see the influence of the evaporation width, we have computed the FWHM 
for ^^''Sn a function of temperature including Ten (evaluated with k = 12), which is shown 
in Fig. 5 by the dashed line. On the whole, the inclusion of Fen leads to a better overall 
agreement with experiment. 

It is to be noted that although the experimental data for ^°^Pb do not, as yet, extend to 
T ~ 3.0 MeV, the effects of the particle evaporation width will also be present in ^°*Pb. We 
have computed Fe,; for ^'^^Pb using the same method as in Ref. [|lOl and is displayed in Fig. 8 
as a function of excitation energy for n = 8, 10, and 12. Also shown in Table IV are values of 
Ybejoreiafter) as a fuuctlou of temperature for k = 10 and 12. The FWHM for ^ospb 

including Fen (with k, = 12) is shown in Fig. 5 with the dashed line, where it is seen that at 
T = 3.25 MeV the FWHM is approximately 2.5 MeV larger than predicted by the adiabatic 
model. As such, experiments carried out in this temperature range would be a further 



confirmation of this effect. Finally, as is pointed out in Ref. |]T0[, the particle evaporation 
width also leads to a maximum excitation energy (or limiting temperature) above which 
the GDR is not observable. This occurs when T^v ~ Fq, which for ^°^Pb corresponds to 
^ 300 - 350 MeV (or T ^ 4.2 - 4.5 MeV). 



VI. CONCLUSIONS 

We conclude that a systematic study of the FWHM of the GDR as a function of tem- 
perature for the nuclei ^^°Sn and ^'^'^Pb confirms, for the first time, the overall theoretical 
picture of the GDR in hot nuclei at low spin. In particular, the role played by adiabatic, 
large-amplitude thermal fiuctuations of the nuclear shape. In fact, overall agreement be- 
tween theory and experiment is observed over a range of temperatures for both ^^°Sn and 
^°^Pb, which display quite different behaviors for the FWHM as a function of temperature. 
This difference can be attributed to the presence of strong shell corrections favoring spheri- 
cal shapes in ^osp]^ that are absent in ^^°Sn. Finally, the increase in the FWHM over that 
expected from thermal averaging at temperatures of the order 3.0 MeV is in accordance 



19 



with the increase expected from the particle evaporation of the compound system. 
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TABLES 



TABLE I. Parameters in Eq. (38) to define the shell corrections to the free energy in ^'^^Pb. 



1 Ai Bi Q 
-13.706 13.764 3.011 

2 -6.448 10.357 3.122 

3 6.68 8.159 2.408 
5 19.50 23.882 3.146 



TABLE II. Parameters in Eq. (41) to define the shell corrections to the moment of inertia in 

208pb. 



I 


Ai 


Bi 


Ci 





-87.653 


13.764 


3.108 


2 


-49.343 


10.357 


3.026 


3 


28.532 


8.159 


3.068 


5 


80.810 


23.791 


3.558 


I 


a« 




Vl 


1 


-22.910 


8.334 


2.979 
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TABLE III. Values of T^n = V^^"'"'' + for ^^OSn obtained from Fig. 2 of Ref. [10] as a 

function of temperature and the level-density parameter defined as a = A/ n MeV~^. All quantities 
are given in MeV 



1 




K = 10 








K = 12 








Ex 


T^fec foTB 

ev 


T~\(iftev 

ev 


Ten 


Ex 


-pijjg f ore 

ev 


T-\ctftev 

ev 


Ten 


1.25 


34 


0.03 


0.00 


0.03 


31 


0.06 


0.02 


0.08 


1.50 


42 


0.06 


0.03 


0.09 


38 


0.10 


0.05 


0.15 


1.75 


52 


0.09 


0.06 


0.17 


46 


0.13 


0.06 


0.19 


2.00 


63 


0.14 


0.08 


0.22 


55 


0.20 


0.11 


0.31 


2.25 


76 


0.22 


0.12 


0.34 


66 


0.33 


0.17 


0.50 


2.50 


90 


0.37 


0.21 


0.58 


78 


0.52 


0.29 


0.81 


2.75 


106 


0.61 


0.38 


0.99 


91 


0.71 


0.49 


1.20 


3.00 


123 


0.83 


0.63 


1.46 


105 


0.95 


0.70 


1.65 


3.25 


142 


1.12 


0.89 


2.01 


121 


1.31 


0.97 


2.28 


3.50 


162 


1.45 


1.19 


2.64 


138 


1.72 


1.36 


3.08 
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TABLE IV. Values of T^n = Tllf"^" + r"!*^"^ for ^ospb as a function of temperature and the 



level-density parameter defined 


as a = 


A/k MeV-\ 


All qunatities are given in 


MeV. 
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"pbe fore 
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Ex 
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^ ev 
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ev 


'- cn 
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61 
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53 
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0.02 


0.09 


1.75 


77 


0.10 


0.04 


0.14 


67 


0.11 


0.06 


0.17 


2.00 


97 


0.24 


0.12 


0.36 


83 


0.23 


0.13 


0.36 


2.25 


119 


0.39 


0.27 


0.66 


102 


0.46 


0.27 


0.73 


z.ou 


1/1/1 


U.D4 


0.47 


1.13 


1 09 


U.DO 


0.50 


1.15 


2.75 


171 


0.90 


0.76 


1.66 


145 


1.01 


0.77 


1.78 


3.00 


201 


1.33 


1.10 


2.43 


170 


1.38 


1.15 


2.53 


3.25 


234 


1.73 


1.52 


3.25 


197 


1.90 


1.58 


3.48 


3.50 


269 


2.20 


2.01 


4.21 


226 


2.42 


2.12 


4.54 
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FIGURES 

FIG. 1. The free energy for ^°^Sn is plotted (lower panel) along the oblate noncollective (/? < 0) 
and prolate collective (/? > 0) axes at a temperature of 2 MeV and a rotational frequency of 
1.25 MeV ((J) fsi 55^). In the upper panel, the Boltzman weight factor exp[— (F — Feg)/r], where 
Fgq is the minimum of the free energy below the saddle point, is plotted. 

FIG. 2. Nilsson-Strutinsky shell corrections (solid squares) to the free energy for ^'^^Pb as a 
function of temperature for oblate (7 = 7r/3), triaxial (7 = 7r/6), and prolate (7 = 0) shapes at zero 
angular momentum. The parameterization to the shell corrections given by Eq. (38) is represented 
by the solid line. 

FIG. 3. Nilsson-Strutinsky shell corrections (solid squares) to the moments of inertia for ^°^Pb 
for 7 = 7r/3, 7r/6, 0, — 27r/3, and — 7r/3. The parameterization to the shell corrections given by 
Eq. (41) is represented by the solid line. The temperature for each panel is the same as in Fig. 2. 

FIG. 4. The weight function W{P) = 13^ /T^l'^e-^''^ at T = 1.0 MeV for prolate (/? > 0) and 
oblate {(3 < 0) shapes. In panel (a), W{fi) includes shell corrections to the free energy, Fshl, 
as well as with (dotted line) and without (solid line) shell corrections to the moments of inertia. 
In panel (b), the same quantities are plotted without shell corrections to the free energy, i.e., 
F = Fld. The free energy with and without shell corrections is plotted in panel (c), and the factor 
{^ld/^SHl)'^^'^ is plotted in panel (d). 

FIG. 5. The FWHM of the GDR strength function as a function of temperature for ^^'^Sn 
and ^^^Pb. Experimental data are represented by the filled circles, while the solid line represents 
the theoretical results obtained for J = 0^. For ^°^Pb, the dotted line is the FWHM obtained 
assuming no shell corrections. For ^^'^Sn and ^'^'^Pb, the dashed line represents the FWHM obtained 
by including the increase to the intrinsic width. Pen, due to the evaporation of particles from the 
compound system. 
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FIG. 6. The FWHM in ^^OSn (dashed hne) and ^ospb (sohd hne) at T = 1.6 MeV as a function 
of angular momentum. 

FIG. 7. The FWHM in ^^Ogn at T = 3.12 MeV as a function of the intrinsic width Fq (sohd 
hne). The experimental value of 11.5 ± 1.0 MeV is represented by the filled square (11.5 MeV) and 
the open circles (±1 MeV). 

FIG. 8. The particle evaporation width, Fev, for ^°^Pb as a function of excitation energy for 
three values of the level-density parameter, i.e. k = 8, 10, and 12 (note a = A/k MeV"-*^). 
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